Influence of disorder on incoherent transport near the Mott transition 
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We calculate the optical and DC conductivity for half-filled disordered Hubbard model near the 
Mott metal-insulator transition. As in the clean case, large metallic resistivity is driven by a strong 
inelastic scattering, and Drude-like peak in the optical conductivity persists even at temperatures 
when the resistivity is well beyond the semiclassical Mott-Ioffe-Regel limit. Local random potential 
does not introduce new charge carriers, but it induces effective local carrier doping and broadens 
the bandwidth. This makes the system more metallic, in agreement with the recent experiments on 
X-ray irradiated charge-transfer salts. 
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I. INTRODUCTION 

Thermodynamic and transport properties of many 
strongly correlated materials are influenced by inhomo- 
geneities and disorder^ and their importance becomes 
even more evident with the advances in experimental 
techniques/?' 3 Very recently, the effects of disorder on 
the optical and DC conductivity of the organic charge- 
transfer salts have been systematically explored by intro- 
ducing defects by X-ray irradiation^— The conductivity 
has proven to be very sensitive on the duration of the 
irradiation, and different physical mechanisms were ad- 
vocated to explain such a behavior4r~ 

Members of K-family of organic charge-transfer salts 
have half-filled conduction bands with the effective 
Coulomb repulsion comparable to the bandwidthji and 
their proximity to the Mott metal-insulator transition 
can be tuned by doping or applying the pressure. On 
the metallic side of the Mott transition, the Fermi liquid 
transport at low temperatures is followed by an inco- 
herent transport at higher temperatures dominated by 
the large scattering rate, and with resistivities an or- 
der of magnitude larger than the Mott-Ioffc-Rcgcl (MIR) 
limit^— which is the maximal resistivity that can be 
reached in a metal according to the Boltzmann semiclas- 
sical theory. From the theoretical point of view, the vi- 
olation of the MIR condition and the appearance of the 
maximum in the resistivity temperature dependence is 
not easy to explain. However, at least for K-organics, a 
significant progress has been recently achieved when the 
transport properties were successfully described even on 
the quantitative level within the dynamical mean field 
theory (DMFT)JJri£ Therefore, it is important to ex- 
plore the effects of X-ray irradiation within the same 
theoretical framework by introducing the disorder into 
the model. 

Most of the theoretical work on the influence of dis- 
order on the physical properties near the Mott tran- 
sition have been so far restricted to binary disorder 
distribution, 16 or low temperatures where the DMFT 



has been extended in order to incorporate the Ander- 
son localization effects*^— However, for a comparison 
with the experiments on K-organics, we need the results 
in a wide temperature range, and since the disorder is 
gradually generated by X-ray irradiation, the simplest 
approach of disorder averaging on the level of coherent- 
potential approximation (CPA)jSi that we apply in our 
work, should be sufficient to explain the main modifica- 
tions in the optical and DC conductivity caused by the 
disorder. We find that, as in the clean case, it is the 
large electron-electron scattering that leads to the de- 
struction of coherent quasiparticles and to the resistivity 
well beyond the MIR limit on the metallic side of the 
Mott transition. For a fixed interaction strength, how- 
ever, the random potential effectively weakens the corre- 
lation effects and moves the system away from the Mott 
transition. Therefore, the disorder increases the conduc- 
tivity in agreement with the experiments on the organic 
charge-transfer salts. 

The rest of the paper is organized as follows. Sec- 
tion II contains the DMFT equations for the disordered 
Hubbard model, brief discussion on the impurity solvers 
that we use and comments on the validity of the DMFT 
approximation. Section III presents the results for the 
temperature dependence of the density of states, opti- 
cal conductivity and DC resistivity near the Mott tran- 
sition for the pure and disordered system. Our results 
are compared with the experiments on X-ray irradiated 
K-organics in Section IV. Section V contains the conclu- 
sions. 



II. DISORDERED HUBBARD MODEL 

We consider half-filled Hubbard model with site- 
diagonal disorder as given by the Hamiltonian 
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Here t{j is the hopping amplitude, U the interaction 



strength, cJ CT is the creation operator, and = c\ a c 
the occupation number operator on site i and for spin 
a. Half-filling condition is enforced by the chemical po- 
tential fi. We model the disorder by random energies 
Vi taken from the uniform distribution in the interval 
(-W/2.W/2). 

In the dynamical mean field theory, which is formally 
exact in the limit of large coordination number, the Hub- 
bard model reduces to a model of an Anderson impurity 
in a conduction bath which has to be determined self- 
consistently. In the presence of disorder, we need to con- 
sider an ensemble of impurities, and the conduction bath 
is determined in the process of averaging over the disor- 
der. 

The central quantity in DMFT is the local Green func- 
tion, Gi G {T-T') = -(Tci CT (r)cJ CT (r'))si , which is a site- 
dependent quantity in the presence of disorder. The local 
effective action is given by 

Slff = -— ^ 4a( ilJ n)[iUn + M ~ v i ~ A(iw„)]c i(T (Uj n ) 
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where A is the conduction bath whose self-consistent 
value will be obtained in the iterative procedure. The 
quantity that we average over the disorder is the local 
Green function, G av (ioj n ) = J dvP(v)G(ioj n ,v). Though 
we consider a continuous distribution of disorder, P(v), 
in practice it is sufficient to take a finite number of ran- 
dom energies, and the integral replaces with a sum. In 
the case of uniform disorder 



f 



N 



G av (iuj„) = Gi(iu} r , 

i=l 



(2.3) 



The averaged Green function G av and the conduction 
bath A determine the self-energy through the relation 

G~v(iu n ) =iuj n + fi- A(iuj n ) - E(iw n ). (2.4) 

The self-consistency condition follows from the assump- 
tion that the lattice self-energy coincides with the impu- 
rity self-energy. Then the disorder averaged local Green 
function has to be equal to the local component of the 
lattice Green function, 



G av (iu) n ) 
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Here D{e) is the density of states in the absence of dis- 
order and interaction. Eq. (|2.4[) determines new conduc- 
tion bath which completes the self-consistency loop. We 
note that our treatment of disorder reduces to the CPA 
approximation in the absence of interaction. 

The most difficult step in the solution of DMFT equa- 
tion is the calculation of the local Green function from 



the Anderson impurity action Eq. (|2.2[) . In this paper, we 
solve the Anderson impurity model with One-Crossing 
Approximation (OCA) , 21 ' 22 and cross-check the results 
with Continuous Time Quantum Monte Carlo (CTQMC) 
impurity solver i 23 ' 24 The OCA impurity solver has an ad- 
vantage that it gives a solution on real frequency axis, 
which is necessary for a calculation of the response func- 
tions. Except at the lowest temperatures, where the 
OCA impurity solver does not reproduce the Fermi-liquid 
behavior, the results obtained with these two impurity 
solvers are qualitatively the same and quantitatively very 
similar. 

Since we are going to compare the results with the 
experiments on K-organics, we briefly comment on the 
validity of DMFT approximation. The main advantage 
of this method is that it treats on an equal footing low 
and high energy part of the spectrum and fully takes into 
account the inelastic electron-electron scattering. There- 
fore, it successfully describes a crossover from the low 
temperature Fermi liquid regime to the incoherent trans- 
port at higher temperatures, which is typical for many 
strongly correlated systems. The DMFT is based on the 
assumption of locality of the self-energy, and hence it 
does not fully take into account intersite correlations. 
However, the K-organics have weakly anisotropic trian- 
gular lattice structure^ which makes the intersite cor- 
relations smaller due to the geometrical frustration. In 
this case, local DMFT approximation gives similar results 
as generalized cluster DMFT.— Finally, since the lattice 
structure enters the DMFT equations only through the 
density of states, the transport properties does not de- 
pend much on the details of the band structure, and we 
will consider the hypercubic lattice which has the density 

of states in the form of a Gaussian D(e) = e~ 2e , 
where the energy is given in units of the half-bandwidth. 



III. CONDUCTIVITY NEAR THE MOTT 
TRANSITION 

The phase diagram of half-filled Hubbard model in 
DMFT approximation is well known/^ At low tempera- 
tures, an increase of the interaction U leads to the first 
order metal-insulator transition. The line of the first or- 
der phase transition ends in the critical point at temper- 
ature T c and interaction U c . In this work we will con- 
centrate on the values of interaction equal and slightly 
lower than U c . At these values of U, there is a gradual 
crossover from the Fermi liquid to an incoherent metallic 
and, eventually, insulating behavior as the temperature 
is increased. 

The central quantity that we calculate is the optical 
conductivity. In DMFT it is given by 2 - 
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FIG. 1: (Color online) Density of states and optical conductivity as a function of frequency in the clean case for U = 0.94 U c \w=o 
(upper panel) and disordered case, U = 0.94 U c \w=i (lower panel). Different colors correspond to the four distinctive transport 
regimes (see the text). The insets show the temperature dependence of DC resistivity. T, ui and W are given in units of bare 
E F . 



Here p(e, v) = — ilmG(e, v), and G(e, v) = (v + /i — e — 
E(i/)) _1 . We will express the conductivity in units of 
the Mott-Ioffe-Regel limit for minimal metallic conduc- 
tivity. The MIR limit, <J MIR , is the conductivity which is 
reached when the electron mean free path becomes com- 
parable to the lattice spacing, I ~ a. According to the 
semiclassical arguments, the electrons can scatter at most 
on every atom and the conductivity in a metal cannot 
be smaller than <J MIR . For half-filled hypercubic lattice 
(which has Gaussian density of states), the MIR condi- 
tion I = a is equivalent to E F r = 1 , where E F is the bare 
Fermi energy, i.e. half-bandwidth of the noninteracting 
electrons, and r _1 is the scattering rate. Here H is set to 
1. Therefore, the MIR limit is set by a condition 

t-) b = -2ImS(0+) = 1, (3.2) 

where S is the self-energy measured in units of Ep . 

The density of states and optical conductivity for a 
clean system and in a presence of moderate disorder, 
W = 1, are shown in Fig. Q] The disorder effectively 
increases the bandwidth and the critical interaction U c . 
In our case, we find that U c \w=o = 2.2 and U c \w=i = 
2.45. The increase of U c due to disorder is in agreement 
with earlier estimates obtained by iterated perturbation 



theoryi^l The critical temperature T c weakly depends on 
the disorder strength, T c \w=i ~ T c \w=o = 0.04, where 
ks is set to 1. On Fig. Q] we compare the data at the same 
relative value U/U c = 0.94, and for several characteristic 
temperatures. We see that the disorder does not lead to 
qualitative differences and if the interaction is the same 
when scaled with U c , the density of states and the optical 
conductivity are even quantitatively very similar. 

We can identify several regimes of the electron trans- 
port. At low temperature (green dotted lines and crosses 
in the insets) the scattering rate, r _1 = — 2ImS(0 + ), is 
small and the transport is dominated by long-lived coher- 
ently propagating quasiparticles. The blue dash-dotted 
lines (blue circles in the insets) correspond to the tem- 
perature when the resistivity is already larger then the 
MIR limit, the scattering rate r _1 is larger than Ep, 
and the Fermi liquid picture of well-defined quasiparti- 
cles ceases to be valid. However, a Drude-like peak in 
the optical conductivity, as well as a peak in the den- 
sity of states, are still present. Our results show that 
this is the case also in the presence of moderate disor- 
der. The resistivity maximum (red full line and square) 
is reached when the peak at the Fermi level is fully sup- 
pressed and when a dip at the Fermi level appears both in 
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FIG. 2: (Color online) (a) Scattering rate (full line) and DC 
resistivity (dashed) as a function of temperature, (b) Quasi- 
particle weight as a function of temperature. The data are 
for the clean system at U — 0.95 U c - 



FIG. 3: (Color online) Temperature dependence of DC resis- 
tivity for different interaction U in the clean case, W = (a) 
and disordered case, W = 1 (b). U is given in units of U C (W). 



the density of states and in the optical conductivity. The 
resistivity maximum is more than an order of magnitude 
larger than p MIR = &~ IR - At even higher temperatures 
(violet dashed line and triangle) low frequency optical 
conductivity increases due to the thermal excitations. 

Fig. [2] helps us to further distinguish the mechanism 
leading to the large resistivity and its strong tempera- 
ture dependence. We see that the scattering rate gives 
the main contribution to the resistivity temperature de- 
pendence and causes the violation of the MIR limit, 
Fig. Ufa), while the quasiparticle (Drude) weight Z = 
(1 + |9ReS(w)/9w| tlJ= o) _1 is almost temperature inde- 
pendent, Fig. [2f b) . The dotted part of the line is an 
extrapolation of the OCA results to zero temperature. 
We have also checked that Z depends very weakly on 
the temperature using numerically exact CTQMC impu- 
rity solver. Therefore, we can conclude that the driv- 
ing mechanism for large resistivity is the large scattering 
rate and not the reduction of the spectral weight near 
the Fermi level. This feature, already seen in the exper- 
iments on VO2 — and charge-transfer salts,— seem to be 
common for the systems with half-filled conduction band 
near the Mott transition. This should be contrasted with 
the doped Mott insulators where the main reason for the 
violation of the MIR condition is a decimation of the 
Drude peak in the optical conductivity by the time MIR 
limit is reached, which can be interpreted as a reduction 
of the number of charge carriers 

The results for temperature dependence of DC resistiv- 



ity, p DC = <j~ 1 {<jj — > 0), for several values of interaction 
U are shown in Fig. [3] The resistivity is given in units 
Pmir- For clarity it is shown on a logarithmic scale. The 
resistivity in the clean and disordered case are even quan- 
titatively very similar when the interaction is scaled with 
U C (W). 

IV. INCREASE OF METALLICITY BY 
DISORDER 

Very recent experiments^ — on the charge-transfer or- 
ganic salts provide a rather unique opportunity to study 
the effects of disorder on transport properties without 
changing external parameters or chemical composition. 
The level of defects (disorder) directly depends on the 
time of exposure to the X-rays. The optical and DC con- 
ductivity are proven to be very sensitive on irradiation 
time showing an increase in the conductivity with the 
time of irradiation. The experiments measured both in- 
terlayer and in-plane resistivity with similar conclusions. 
Different physical mechanisms were proposed to explain 
the increase of conductivity. Analytis et al.^ proposed a 
defect-assisted interlayer conduction channel for the re- 
duction of resistivity, and Sasaki et al^^ proposed that 
the irradiation leads to the effective doping of carriers 
into the half-filled Mott insulator. 

The DMFT has successfully described the transport 
properties of organic salts even on the quantitative 
leveL 14 ' 15 In order to make a comparison with the ex- 
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FIG. 4: (Color online) (a) Temperature dependence of DC 
resistivity for fixed U = 2.2 = U c \w=o, and various levels 
of disorder, (b) Experiments on k-(BEDT-TTF) 2 Cu(SCN) 2 , 
taken from Ref. U 



perimcnts with irradiation induced defects, we solve the 
DMFT equations for fixed interaction U and vary the 
level of disorder W. The results for DC resistivity are 
shown in Fig. 0Ja). The data for T < 0.01 are obtained 
using CTQMC impurity solver. The presence of even 
a weak disorder significantly decreases the resistivity by 
effectively moving the system away from the Mott insu- 
lator, as explained in the previous section. Our data are 
very similar to the measurements on charge-transfer salt 
k-(BEDT-TTF) 2 Cu(SCN) 2 from Ref.!, which are shown 
in Fig. Bib). We note that these data are for interlayer 
resistivity while our DMFT calculation corresponds to 
in-plane transport. However, the interlayer transport is 
due to incoherent tunneling which is proportional to in- 
plane conductivity^ Therefore the temperature depen- 
dence of out-of-plane resistivity should follow the tem- 
perature dependence of in-plane resistivity. Indeed, the 
in-plane optical conductivity measurements on the Mott 
insulator rc-(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl, also show that 
the Mott system becomes more metallic in a presence of 
disorder. These measurements show the transfer of the 
spectral weight to low frequency region as the irradia- 
tion time increases, followed by the collapse of the Mott 
gapj£i& 

We emphasize that our model, as opposed to the phys- 
ical mechanism proposed in Ref. 6, does not assume an 
introduction of new charge carriers since the total num- 
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FIG. 5: (Color online) Local occupation number per spin as a 
function of site disorder m for T — 0.01, U = 2.1, and W = 1. 



ber of carriers per site remains equal to one. The lo- 
cal occupation number, however, depends on the ran- 
dom site potential, and we can say that the system is 
effectively locally doped<2£ The occupation number, for 
a given spin orientation, as a function of random site po- 
tential is shown on Fig. [5] It is interesting to note that the 
local occupation number, n(Vi), deviates much less from 
its average value than it would be the case in the absence 
of interaction. This is a consequence of very strong disor- 
der screening of site-diagonal disorder on the metallic side 
of the Mott transition.— Therefore, the resistivity curves 
on Fig. EJa) cross at very low temperatures and our cur- 
rent model cannot explain the intersection of curves in 
Fig. 0Jb) which happens at much higher temperature. 
The dramatic reduction of the elastic scattering is also 
demonstrated in Ref. |32|, which shows that the inelas- 
tic scattering dominates in the incoherent regime. We 
stress that we do not assume Matthiessen's rule. This is 
a salient feature of DMFT, which can operate in a regime 
where conventional approaches to the electron transport 
fail. 



V. CONCLUSIONS 

In summary, we have examined the influence of random 
potential on the optical and DC conductivity for half- 
filled Hubbard model in a vicinity of the Mott transition. 
Our results show, in agreement with the experiments on 
K-organics, that the disorder can make the system effec- 
tively more metallic. The disorder increases the band- 
width which increases U c and weakens the correlation 
effects, moves the system away from the Mott transition 
and leads to a decrease in the scattering rate and resis- 
tivity. We emphasize that the randomness in our model 
does not change global doping, as the system remains on 
average half-filled, but the number of charge carriers lo- 
cally deviates from the average value. Therefore, global 
carrier doping of a Mott insulator due to irradiation de- 
fects, proposed in Ref. [|1 is not necessary to make the 
system more metallic. We also find that the maximal 
possible value of metallic resistivity remains more than 
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an order of magnitude larger than the MIR limit even in 
a presence of moderate disorder. As in the clean case, 
the violation of the MIR limit is driven by a large scat- 
tering rate due to the electron-electron scattering, and 
Drude-like peak in the optical conductivity persists even 
at temperatures when the resistivity is well beyond the 
MIR limit. 
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